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In this paper analysis of two radiation problems which have practical application to improve radio 
communications during re-entry blackout period is presented. In both of these examples given 
here, the mutual orientation of the antenna and an external static magnetic field is chosen in such a 
way that the field components are independent of the component of the plasma-dielectric tensor parallel 
to the static magnetic field. This choice enables one to control more effectively the electromagnetic 
waves by controlling the applied d-c magnetic field. 

Though the problems investigated here involved cylindrical body of infinite length, a knowledge 
of these results will throw more insight into the expected behavior of the radiation field from a uniform 
magnetic ring current around a conical space vehicle covered by a plasma sheath in the presence of 
a uniform azimuthal static magnetic field. 



1. Introduction 

Radio communication with all types of manned space vehicles is desirable during all phases 
of the flight mission. When a space vehicle enters the earth's atmosphere at hypersonic speed, 
a plasma sheath is formed around the vehicle and the antenna installed on it. The electron con- 
centrations in the plasma sheath thus formed may be high enough to degrade communications with 
the vehicle. 

There is also the familiar phenomenon in connection with radio wave propagation through the 
ionosphere that in highly ionized media (which is characterized by the plasma-electron angular 
frequency o) p ) a comparatively low frequency signal (i.e., cu < a> p , oj being the radian signal fre- 
quency) cannot propagate. Use of high frequency (co > cu p ), on the otherhand, does not always 
seem to be very practical, because the available power at such a high frequency is not sufficient. 
However, it is also a known fact that radio wave propagation in the ionosphere again becomes 
possible (even when w < a> p ) whenever it is properly influenced by the earth's static magnetic 
field. Therefore, it is very natural for one to borrow this idea of somehow introducing a uniform 
static magnetic field in the plasma sheath surrounding the space vehicle during its re-entry phase. 
Actually, the idea of using a static magnetic field is not new, and there is a good deal of discussion 
on this subject in literature [Papa and Allis, 1961, and Hodara, 1961]. 

All of these discussions, however, pertain to simple geometry and analysis of plane wave prop- 
agation in magnetoplasma media. Therefore, effects of boundary of a space vehicle on the prop- 
agation of electromagnetic waves in a plasma were not taken into account. Unfortunately, most 
of the space vehicles of practical interest have complicated shapes and an exact analysis of elec- 
tromagnetic wave propagation affected by such boundaries is beyond the scope of the present state- 
of-the-art. However, it is worthwhile to investigate the possibility of re-entry communication with 
a space vehicle having simple geometries, such as a semi-infinite conducting cone or an infinite 
conducting cylinder. 
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Using such simple configurations, it was observed qualitatively in a recent communication 
[Samaddar, 1965] that a satisfactory elimination of the blackout depended considerably also on the 
mutual orientation of the static magnetic field in the plasma sheath and the antenna. In particular 
it was proposed that the orientation of the antenna and the static magnetic field should be chosen 
in such a way that the field components were independent of the component of the dielectric tensor 
parallel to the static magnetic field. This choice enables one to control more effectively the propa- 
gation characteristics by controlling the applied d-c magnetic field alone. 

In this paper analytical results of the above mentioned proposed problems involving cylindrical 
geometry are presented. Since the special functions used in the analysis are not well tabulated, 
the numerical results are left to another paper in the near future. 

The first problem considers the radiation from a magnetic current line source (an idealization 
of an axially slotted antenna) on an infinitely long metal cylinder, covered by a cold plasma in an 
axial static uniform magnetic field. In the second problem the study of wave propagation through 
a plasma-clad metal cylinder in a uniform, angular, static magnetic field is made. The source of 
excitation is a uniform magnetic ring current situated coaxially on the surface of the metal cylinder 
so that it will excite only the axially symmetric E-type modes. 

It may be mentioned here that in both of the above mentioned problems, the electromagnetic 
fields do not depend on the component of the dielectric tensor (for the plasma) parallel to the static 
magnetic field [Samaddar, 1965]. Note that if this mutual orientation of the antenna and the static 
magnetic field is not preserved, a complete elimination of the blackout by simply increasing the 
strength of the static magnetic field cannot be attained in principle. 

In passing it may be pointed out that though a theoretical study of the radiation of electromag- 
netic waves from a uniform magnetic ring current around a conical space vehicle covered by a 
plasma sheath in the presence of a uniform angular static magnetic field is very difficult, the results 
of the second problem mentioned in the preceding paragraphs will throw more light into the prac- 
tical aspect of the conical configuration. 

2. Problem Involving Axial Static Magnetic Field and Magnetic Current 

Line Source 

2.1. Statement and the Formulation of the Problem 

Consider a space vehicle which can be represented, for certain practical purposes, by a con- 
ducting infinite cylinder of radius a, which is surrounded uniformly by a plasma sheath of radius 
b in presence of a uniform axial static (d-c) magnetic field (see fig. 1). A magnetic current line 
source is situated axially along the surface of the metal cylinder. The region outside the plasma 
sheath is assumed to be free space. Due to the cylindrical symmetry of the problem, the cylindri- 
cal coordinate system p, 0, and z will be used. Thus, the coordinates of the source are given by 
a and Oq. 

The method of procedure to solve this problem is already known [Ohba, 1963; Samaddar, 
1962; Wait, 1961]. As a matter of fact, the formal result is already given by Ohba [1963]. There- 
fore, we shall merely present the result here on the basis of the procedure adopted by Samaddar 
[1962]. In addition to the formal result, a few special cases will be derived using appropriate 
limits. 

Since the cylindrical structure is independent of the z-coordinate (axial direction) and the 
source is a uniform magnetic current line source parallel to the z-axis, the field components (Ep, 

Pi 

E$, and H z ) will be independent of the z-coordinate, i.e., — = 0. In a straightforward manner, 

oz 

it can be shown that the longitudinal field H z (assuming a suppressed harmonic time dependence 
e t(ot ) has the following representations: 1 



1 Note that the field components in this radiation problem are independent of e zz , the component of the dielectric tensor parallel to the d-c magnetic field B . 
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Figure 1. A magnetic current line source on the surface of an infinitely long metal cylinder covered by a magnetoplasma 

sheath. 



Hi • ^Sf ? ~^~ Kki^WW-tf. 6 > P)-L«Qn({; b,p)], for a ^ p =£ b (l a ) 



where, 



ima)e ei6' ^, e im ~ e o> 



Hf(K p)iorp>b 



W„ = L n R n {t- a, b) + #>e,S»(£; a, b)H™(K b), 
e' = [e?-e*]/e, 

(1)1(1 — iv/(o) 



€l — € pp — 600—1 + 



(o 2 — (o 2 (l — iv/co) 2 



(lb) 

(2) 
(3i) 

(3ii) 

(3iii) 

(3iv) 



€p0 — 6o P — €2 : 



(O c 



co (o\— co 2 (l— iv/co) 2 



top = angular electron-plasma frequency = a/ » 



e = magnitude of the charge of the electron 
No = number density of the electrons 
m e = mass of an electron 



e = free space permittivity 
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(3v) 

(3vi) 

(3vii) 

(3viii) 

(3ix) 

(3x) 



(i) c = = electron cyclotron frequency (angular) (3 x j) 

B = static (d-c) magnetic induction in the z-direction (3xii) 

m = strength of the magnetic current per unit length (3xiii) 

v = collision frequency (3xiv) 

L n = £bV^ e 1 H^\Kob)-n€ 2 H^(Kob) ( 4i ) 

HWKob)=-^H%KKob) (4ii) 

Rn(i; a, b) = J n (tb)N'Ma)-J'Ma)N n m (4iii) 

S n ({; a, b) = J' n {ta)N' n {Zb)-J' n (Zb)N'Ma) (4iv) 

Q n (t; 6, p)=M£p)Nn(€b)-Mtb)Nn(£p) (4v) 

R n ({; 6, p)=Jn{i;p)N t n m-J' n mN n {Zp). ( 4v i) 

2.2. Special Cases 

In absence of the static magnetic field, the plasma becomes isotropic. In this case, the ex- 
pression for the field component H z can be calculated from (1) by employing the following limits; 

#o = 0, *2 = 
and 

/ -. <°l 

e'=€ 1 = € p = l- v-. • (5) 



Using these limits of (5), the expressions for the field component H z can be written in the 
following form: 

ft = imo)6o6 p J 5__ [66ff (2) (Xo6)/?n(6;6>p) 

-gtbV^HWiKobWnitiib-p)] 

for a< p < b (6a) 



and 



where 



rr im G)€ e2 • c *n(0-0o) 
^ n=-oo U1 



(6b) 



r ln = ^be P [V7 p H^\K b)RM, ;a,b)+ S»(f , ; a, b)H™(K b)]. (?) 
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When there is no plasma around the conducting cylinder, the expression for the field can be 
obtained from (6) using the limit a = b and e p = l. Then one finds 



= imve l) //i, 1) * _j 
lira .A H 



..;r(K a) 



H<*\Kop), p>a. 



(8) 2 



For a thin plasma, i.e., for (b — a) < a, £(b — a) < 1 and K a (b — a) <§ 1 [Wait, 1959], the expres- 
sions of H z for p > b can be obtained by Taylor's series expansion of the cylindrical functions and 
their derivatives in the following manner: 



and 



Z n (ria) = Z n ( V b) + 7)(a - b)Z' n (r)b) 



Z' n (r,a) = Z'^b) + V (a-b)Z:(r,b) 



(9a) 



(9b) 



where Z n is any of the cylindrical functions and r) may represent £, fi, or K as the case may be. 
In doing so, it is a simple matter to show that 



R n (y; a, b) 



TTTja 



(10a) 



W n 



TT 



S n (r}; a, b) ^ 



2(6 -q) A n 2 



nb 



V V 2 bV 



- a V^e l H'J»(Kob)+ { 6 »« & -«)( 1 -^)-^ L } U*\K b) 



(10b) 



(11a) 



and 



Wu 



2VT V > 



'& [- a H'*\Kob) + Xo(6 - a) (l - ^) flW^W)]- 



(lib) 



Now substituting (11a) and (lib) into (lb) and (6b) respectively, we obtain the corresponding 
fields for the thin plasma sheath approximations. 

From a knowledge of the numerical results obtained previously [Ohba, 1963; Samaddar, 
1963-1964; Seshadri, 1964] for problems of this nature, some general remarks can be made here. 
Both the amplitude and phase of the far field behave asymmetrically with respect to the angle 6. 
This asymmetry is caused by the presence of the static magnetic field in the plasma. However, 
if the excitation is such that the electromagnetic fields do not experience [Samaddar, 1962] the 
anisotropic nature of the plasma, no asymmetry of the amplitude and phase as mentioned above 
will occur. For some given angle 0i, the amplitude may increase and the phase decrease with the 
increase of the static magnetic field. However, this increase of the amplitude of the radiation 
field with the static magnetic field may not be true for all angles 6, \6\^tt. This suggests that 
the application of an external magnetic field may change both the gain and the directivity of the 
radiation pattern. 



2 Note that this result can also be obtained from (1) by taking the limit a> r /u) ► 1 and a> r /a> p > 1. In other words, in this limit, re-entry blackout disappears. 
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Figure 2. A magnetic current ring source on the surface of an infinitely long metal cylinder covered by a magnetoplasma 

sheath. 



3. Problem Involving Angular Static Magnetic Field and Magnetic Current 

Ring Source 

3.1. Statement and the Formulation of the Problem 

In this problem we study the radiation of electromagnetic waves from a uniform magnetic 
current ring source (i.e., an idealization of a circumferential slot antenna) on the surface of an in- 
finitely long conducting cylinder. The cylinder is covered by a plasma sheath concentric with the 
cylinder. An external static (d-c) magnetic field in the angular direction is assumed to exist in the 
plasma sheath (see fig. 2). The space outside the plasma is unbounded free space. 

We further assume that the interaction of the plasma and the electromagnetic fields can be 
described by using the following linear hydrodynamic theory of a cold plasma together with the 
Maxwell's equations (with a harmonic time dependence e i(ot ): 



V X E =- iw/xoH -OoM (12a) 

VxH = i(oe E - eN \ (12b) 

m e (v + ico)\ = — eE + eB 0oX V, where 0o is the unit vector in the angular direction. (12c) 



From (12b) and (12c), it is evident that the motions of ions are neglected in comparison with those 
of the electrons. V is the a-c velocity of the electrons, and all other symbols have 'heir usual 
meaning. The magnetic current ring source M can be represented by Dirac delta functions in the 
following manner: 



M~- 



md(p — a)5(z) 
2?rp 



(13) 



where m is the strength of the source and a is the radius of the cylinder. 

Due to the choice of the source of excitation, the electromagnetic fields will be independent of 

a 



the angular coordinate 6 (i.e., 



dO 



'■ 0). Now in a straightforward manner, it can be shown that only 
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the components E p , E z , and He will be excited and they obey the following partial differential 
equations: 



dE p dE z . rr mS(p — a)8(z) 



— - = — i(o€ (€iE p + ieJEz) 
oz 



l_d_ 
pdp 



(pHe) = i(D€ (eiE z — ieJLp). 



(14a) 
(14b) 
(14c) 



where €i and e 3 are the components of the following dielectric tensor: 



e = 



€i 





l€ 3 





€ 2 





i€z 





6] 



(15) 



€i = €p P = €j W has the same representation as given by (3iv), and €3 = € pz = € Z p has the same 
expression given by (3v), whereas 62 — eee can be shown to equal to 



1- 



(o 2 (l — iv/(o) 



It may be noted here that the field equations (14) representing this radiation problem do not depend 
on the component, 62, of the dielectric tensor. Note that €2 is the component of e parallel to B . 
In order to solve the partial differential equations in (14), it is convenient to take their Fourier 
transform with respect to z, i.e., if P(p, z) is any component of the field, then its Fourier transform 
is defined by, 



P(p,z)=^["P(p,£)e-«*d£ 



(16a) 



P{p, 0=\ P(p, z)e«*dz. 

J -00 



(16b) 



Now employing this operation of Fourier transform in each of the equations in (14), it can be shown 
after some rearrangement that the transformed field components obey the following ordinary 
differential equations: 



d 2 * 1 d A 
-7-5 He H ~t~ tie + 
dp 1 p dp 



€l p 2 e,pj 



= ia)€ m * 8(p — a) (1 7a) 



2jrpe\ 



E z = 



1 



KO€o(€f — €§, 



r ia 

L ei p dp 



{pH e )-ieMe 



(17b) 
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It may be noted here that in absence of any static magnetic field, e^ = 0, then the homogeneous 
differential equation corresponding to (17a) could have been satisfied by Bessel and Hankel func- 
tions. However, in this present problem, (17a) can be satisfied by two linearly independent 
confluent hypergeometric functions and the required solution for He can be written in the fol- 
lowing form: 



H*p, Q=A (af ) r*» (^±£ 3, 2jqp) 

-B ^L (top)**** (^±£, 3, 2fr|p) 

= CH?X. Vo p), for p&b, 
where A, B, and C are arbitrary constants and 

V = ^(ef-e2)/ €l -C 2 
f-fe/ei. 



for a =£ p « 6 



(18a) 



(18b) 



(19) 



this £ should not be confused with that used in section 2. 

The confluent hypergeometric functions <l> and ^ used here are in accord with the definitions 
given in Bateman Manuscript Project, [1953] and Slater [I960]. The relation between the singular 
solution ¥ and the regular solution <P can be expressed in the following way: 



¥(a, 3, 2ir,p)=- 



1 



2IXa-2) 



<D(a,3, 2ir,p)log(2i7,p)+yg r 



{«fta + r)- 



where, 



W+D-^}«| +i L 



3, if 

a = 2 + 2^' 



\2i7jpj 



+ 



2ir)p 



d 



My)=^ r(y)/r(y). 

(a) r = r(a + r)/r(a) = a(a+l)(a + 2) . . . (a+J—1) 
(3) r =r(3 + r)/r(3)= 3.4.5 . . . (3+r-l). 



(20) 

(21a) 

(21b) 

(21c) 
(21d) 



The main reason for writing the equation (18a) in a little cumbersome way is to show that in 
the limit of £ = 0, (i.e., in absence of the d-c magnetic field), the functions <£ and ^V reduce to Bessel 
and Hankel functions, respectively, as given below: 



Mv,p) = ij f-e-« 



<D (| 3, 2n,p) 



and 



Hf\ V p)=-^=(2r,p)e-^^ (| 3, 2n,p). 



(22a) 



(22b) 
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Though it looks elegant to use the known confluent hypergeometric functions <f> and ty as the solu- 
tions of (17a), for numerical computations this representation may not be convenient, since the 
functions <J> and ^ with complex arguments are not tabulated. Therefore, it may be worthwhile 
to express <J> in terms of the Bessel function J\(j)p) and a series as given by (23). Once <J> is calcu- 
lated from (23), the singular solution *P can be computed by using (20) and (23): 



*(l + V 3 ' 2 ^) = 4 e, '" p [Mvp)+id n (vp) n+t ], 



(23) 



where 



and for n ^ 2 



<fe = f(ll7f-|*)/(48.15ij s ) 



(24a) 
(24b) 
(24c) 



4-ii){(2n+ l) 2 - 1} + d 2 „- 2 f + di„. :! rj=0 

(-D" +i g 



«fe»T?{4<n + D 2 - 1 } + &»- if +■ A„- 2 17 : 



2 2n+l ra!(re+l)! 



(24d) 
(24e) 



Note that d n = for f = 0. Before calculating the unknown coefficients .4, B. and C by using the 
appropriate boundary conditions, we present the following relations which will be found very 
useful in the sequel: 



The Wronskian, <t> — V-ty — (J> = - 2(2njy)~ *««w/r(o) 



1 d 



*'W = 2T^-*(«. 3 > 2JT)y)=-a¥(a + 1, 4, 2njy), 



(25a) 
(25b) 



*'(y)=7T"T-*(a, 3 > 2»Tiy)=^*(a+l, 4, 2»wy). 

Zit) ay o 



(25c) 



The transformed field //»(p, £) given by (18) satisfies the following boundary conditions: 



e?-«S) 



He(b-,Q = H e (b\i\ 



1 J « * 1 1 w 

ei~i fi H t )-i€ t H l 

pdp 



i<*™ 



p=ft - p dp 

_ imo)6o(e 2 — e 2 ) 

2776! 



ipHe) 



p=6" f 



(26a) 
(26b) 

(26c) 



The relation (26c) specifies the source condition. Now applying the boundary conditions 
(26) on the solution (18) together with the relations (25), the transformed field component Heip, Q 
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can be shown to be equivalent to the following expressions: 

*<* 0= ^^ N( r(OT) -^}- r(a)W H < 27a > 

for a ^ p ^ 6 

HslP ' g 4 7 r(7 ) a)(T ? 6)e 1 7/( 2 >(T,o6)[(?i - <M ' l W ' P (27b) 

<D(p) = *(a, 3, 2iTjp), ¥(p) = ¥(a, 3, 2i7jp), 

<?, = [(2 - inaMa) + 2h)a&(a)] \T(a)V(b)Q -^jTJ • (28a) 

<? 2 = [(2 - i V ama) + 2i>qaV'(a)yS>W(a)Q. (28b) 

^K fi ir l ) ^S " (2 ~ ** ~ ir,b) ] ~ 2 W<*>> (28c) 



and 



where, 



and 



Q = <P(b) 



^'(a) and 0'(6) are defined by (25b) and (25c) with y = a or b as the case may be. 

Since we shall be interested in the far field, we shall now take the inverse transform of He{p, Q 
only for p> b. Therefore, taking the inverse transform of (27b), we have 

-imc6o(6?-e 2 ) f- e^^r,^-^, 

»»</>> *J 87r 2 €ia6 J-^ * v'HfKvMQi-Q*] (29) 



3.2. Evaluation of the Integral lor the Radiated Field 

Now defining a new function F(Q by 

-ima>eo(e 2 -e 2 )<&(6)e'""' + '-> 
r ^ fb*€ ia bri*H<*K7,MQi-Q2] ( ' 

the integral in (29) can be written in the following way: 

H<tp, z) = | y F(OH^( V op)e-^dC (31) 

where the path of integration y runs along the real axis of the complex £ plane with identations at 
the branch points £=±£ and at the real poles (if any) in such a way that the far field satisfies 

1/2 

the radiation condition at r = (p 2 + z 2 )^ oo. In other words, the integral converges (fig. 3). 
The poles of F(£), which corresponds to surface waves are given by 

Qi = Q2. (32) 

Since these surface waves, unless they are encountered with any discontinuity along the direction 
of propagation, have no appreciable influence on the radiated field, they will not be considered 
here. A detailed study of the surface waves for this problem is left to a subsequent paper. 
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FIGURE 3. Path of integration, y, and branch cuts in the complex £ plane. 



Though 7) which appears in the integrand, is a multiple-valued function of £ in the neighbor- 

/ € 2_ € 2\ 1/2 

hood of £ = ± K I — ) it can be shown that the integrand is an even function 3 of 17. Therefore, 

/ e 2_ e 2\ 1/2 

17 = 0, i.e., £ > = ±Ko(— -) are not branch points of the integrand. However, at 170 = 0, the 

Hankel functions have logarithmic singularity. Therefore, £=±X (i.e., 170 = 0) are branch points 
of the integrand. The figure 3 shows the choice of the branch cuts. 

Once the singularities of the integrand are identified, the evaluation of the integral (31) by 
the steepest descent method for rjop ^ 1 is a straightforward procedure (assuming that there is 
no pole near the saddle point, nor on the steepest descent path). Therefore, without going into 
detail the first order asymptotic value can be given by 



H e (p,z) 2F(£ = X cos/3) 



-, for T7op > 1. 



(33) 



where 



z — r cos ft 
p — r sin f3 
C = K cos /3 
170 = K sin /3 

and r, /3, 6 define a spherical coordinate system. The function F(Q of £ is replaced by F{K Q cos (3). 



3.3. A Few Special Cases 

In absence of any static magnetic field in the angular direction, B = 0, e.3 = 0, £ = 0, ei = e p 
= l — (D 2 l[(o 2 (l — ivla))], and using these limits together with the identities (22), the following rela- 



3 To show this, it is useful to know that e' ivp 4>( \ % 3, 2ir)pjis also an even function of tj, which follows from (23). Furthermore, the behavior of rjp e ' 



^l ^9 ' ^' 2m?P )as tj approaches zero is similar to that of H\ v (r)p)- 
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tions can be established: 



B» Q> = inei Zw(v!t a) ^beM\vob)Hf\ v b)- v bH^(vb)HfKvob)] 



Bo-+0 



lim Q 2 = ^ e ' 2(^ff ) >(T?a) l-nobepJi(vb)H*\r to b)-'nbMr,b)H?\r,ob)] 



B ^ 






flo-*0 



L, = r,bHM(Vob) [Mva)No(-nb) - /ofatyVodja)] . 
L 2 = r, be p Hf(y) Jb) [MwWiiyb) - y.faWVofaa)] . 



(34a) 
(34b) 

(34c) 

(34d) 
(34e) 



with 



Now, using these limiting results, we have 

lim FtO-Fi®- *""** . V 

s„-*o 2TT J ria(L,—L2) 

Therefore, the far field expression in this case becomes 

e -iK r 

Hdp, z, Bo = 0) = H l0 {p, z) ~ - 2 Ft(i = K cos j8) — — , for ij p > 1. 

In the limit of a thin plasma sheath, i.e., T)(/> — a) <§ 1, with Bo — 0, Fi(£) becomes, 

imcuepep 



Fi(Q- 



ktth\a Ui(b - a)Hf>( V ob)+^e p H%\ Vo b) 



(35) 



(36) 



(37) 



For the situation where there is no plasma sheath surrounding the conducting cylinder, one finds: 

a = b, € p = 1 , Z-i = 0, 17 = 170 



U = --HU\r) a) 



lmcueo 



.-. lim F(Q = lim F,(0 = F 2 (£) = - 

*-o .,-1 47r 2 i7oa//< 2 )(T7oa) 



(38) 



Therefore, the far field is given by 



Hdfi, z, B - 0, e„ - 1) - ff 29 (p, z) 2 F 2 (£ = Ko cos /3) - 



(39) 



In the hmit of a very thin plasma sheath, FiiQ can be expressed in terms of F 2 (£) in the follow- 
ing way: 



Fi(Q~FdQ 



1+ u?2W — TsT~ ' for 7)(6 - a) « 1, 



(40) 
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where A = 1 — ifteplrf-. 

It may be noted here that if one uses the limits co r /o> > 1 a> c /a> p !> 1 carefully in (27) and (33), 
one can show that (33) reduces to (39). This implies that (even without performing numerical 
computation) in the limit of high d-c magnetic flux, the plasma behaves as the free space and thus 
the blackout hazard disappears. 
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